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We present the first high-redshift gal axy cluster candidate sample from the HIROCS survey found in the COSMOS field. It results 
from a combination of public COSMOS with proprietary H-band data on a 0.66 part of the COSMOS field and comprises 12 
candidates in the redshift range 1.23 < z < 1.55. We find an increasing fraction of blue cluster members with increasing redshift. 
Many of the blue and even some of the reddest member galaxies exhibit disturbed morphologies as well as signs of interaction. 
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Introduction 



Galaxy clusters provide ideal laboratories to study galaxy evo- 
lution. Furthermore, they can also be used to constrain cosmo- 
logical parameters ( |Voit||2005] l as well as theoretical models 
of environmental effects for which there is growing evidence 
( |Khochfar & Ostriker|2007| . Only few galaxy clusters at z > 1 
are known today, most of them stemming from X-ray surveys 
(Mull is et al.||2005[ [Stanford et al.|2"0"06*l ). These detect clusters 
with a hot intra-cluster medium, as will future surveys based on 
the Sunyaev-Zel'dovich effect ( |Kneissl et al.|200T] >. Optical se- 
lection of clusters provides a complementary sample ([GTlbank 
et al.||2004 i, tracing stellar light. At z ~ 1.1, where the 4000 A 
break gets shifted to the near-infrared, optical surveys need to 
be supplemented with at least one near-IR band. With the advent 
of wide-field near-IR cameras such as OMEGA2000 on Calar 
Alto, wide-area imaging surveys of z > 1 clusters have become 
feasible (Elston et al. 2 006{ [Lawrence et al.|2006| . A number of 
optical / near-IR surveys is now beginning to find clusters in this 
redshift range (Stanford et al. 2005] [van Breukelen et al.|20 06 ). 

In this letter we use proprietary //-band data in combina- 
tion with COSMOS public data ( |Scoville et aLl|2"0"0"6[ > to find 
clusters in the redshift range 1.2 < z < 1.6 in the COSMOS 
field. Our detection method employs photometric redshifts to 
estimate the distances and finds clusters as overdensities in 3D 
galaxy distribution. Hence, not all our candidates need to be viri- 
alized. Some of them may also be clusters in the state of forma- 
tion. Throughout this letter we assume H = TOkms-'Mpc- 1 , 
£2 m = 0.3 and Q\ = 0.7. All magnitudes quoted are in the Vega 
system. 



Send offprint requests to: M. Zatloukal, e-mail: Zatloukal@mpia . de 
* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck-Institut fur Astronomie and the Instituto de Astrofisica de 
Andalucfa (CSIC) 



2. Data 

HIROCS, the Heidelberg InfraRed / Optical Cluster Survey, 
aims at finding galaxy clusters at 0.5 < z < 1 .5 by means of pho- 
tometric redshifts. The survey area is 1 1 n° in four fields, one of 
them including the COSMOS field. A detailed description of the 
survey will be found in Roser et al. (in preparation). 

In the COSMOS field, HIRO CS uses the public data in u* ', 
Bj, g+, Vj, r+, i+, NB816 and z+ ( |Capak et al.]2007] >. K s is used 
only as a supplement since it is very shallow compared to the 
other bands. We use the spatially psf homogenized matched_psf 
mosaics with PSFs from 0.9" to 1.6" for the individual fil- 
ters. The data come calibrated in units of nJy with informa- 
tion about errors and bad pixels stored in separate RMS images. 
Morphological information is provided by HST ACS F814W 
data ([Koekemoer et al. 2007]l. 

As a complement to the optical COSMOS data we use pro- 
prietary //-band data taken with OMEGA2000 at Calar Alto 
Observatory in a joint effort of the large extragalactic sur- 
vey projects ALHAMBRA ( |Moles et al.|2005| l, MUNICS-deep 
(Goranova et al., in preparation) and HIROCS. The COSMOS 
field will be covered by 25 OMEGA2000 pointings of 15.4' x 
15.4' with a targeted depth of 21 mag (5<x) and seeing < 1.6". 
Currently, the //-band data cover 1 .3 with a non-uniform 
depth. For this letter we use the deepest 0.66 n° which reach 
21.4 mag (3cr) or deeper, see Fig. [2] 

3. Data reduction and analysis 



The OMEGA2000 reduction pipeline ( |FaBbender| [2007 ) was 
used to reduce the //-band data applying object masking for 
sky determination and image stacking using relative weights and 
fractional pixel offsets. Flatfields are derived from the science 
data. Bad pixel and cosmics with a difference of > 3<x to a me- 
dian image are replaced by the scaled median value. The images 
are astrometrically aligned using point sources from the 2MASS 
catalogue (Skrutski e et al.|2006| >. 

Our final object catalog contains all objects from all the 
bands regardless of their colour with positions as input for the 
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photometry. We use SExtractor (Bertin & Arnouts 1996 ) to ex- 
tract an object list from the coadded images for each filter with 
a minimum S/N threshold of 5, omitting the shallow K s images. 
Spurious detections, e.g. in the halos around bright stars, are re- 
moved manually. Next, an artificial image containing all objects 
from the different SExtractor tables is created. For each object, 
a Gaussian is added into this image with the S/N as amplitude 
and the shape as measured. The final catalog is created by run- 
ning SExtractor on this artificial image with a special setup for 
noise-free data. It reaches a mean positional accuracy of ~ 0.1". 

We use MPIAPHOT to derive instrumental magnitudes as 
weighted sums over the image area (Roser & Meisenheimer 
|1991[ l. The width of the weighting function is adjusted such 
that the measurements refer to a common Gaussian "aperture" 
of 1.8" for all images: cr^ nmmnn = <x 2 . . , + a 2 ■ , equivalent 

to common weight seeing' " 

to the beam in aperture synthesis radio astronomy. Photometry 
in the COSMOS bands is done on the coadded frames. The flux 
errors are measured using the same parameter files on variance 
frames created by squaring the flux of the RMS frames, //-band 
photometry is done on subsums of 5 one-minute exposures coad- 
ded using relative weights and fractional pixel offsets. 

The //-band needs a mosaic correction of the instrumen- 
tal magnitudes due to varying observing conditions. Each of 
the 25 pointings is calibrated using 2MASS stars brighter than 
H = 14.5. We get a mean zeropoint error of 0.008 with an 
RMS of 0.07 measured across the whole field. Fine tuning of the 
photometric calibration zeropoint, especially to relate the optical 
COSMOS magnitudes to the infrared magnitudes, is done by ap- 
plying a stellar main sequence colour shifting technique similar 
to the one used by | Wolf et al^pOOT) . 

Our object catalog comprises 32798 objects with H < 21.4 
to be subject to multi-colour classification. Due to the deep 
COSMOS optical data, all objects have photometric redshifts as- 
signed and our catalog is essentially complete down to H — 21 .4. 

The multi-colour classification and photometric redshift es- 
timation is done using the template-based code developed for 
CADIS and COMBO- 17 by [Wolf et alT| < |2001> . The quality of 
the photometric redshifts is demonstrated in Fig.fTI We use 140 
spectroscopic redshifts from Pres cott et al.| (|2006| leavin g out 
the QSOs and a sample of 48 zCOSMOS ( |Lilly et al.|2006) ob- 
jects from masks obtained from the ESO archive (Gabasch et al., 
in preparation) to compare our photometric redshifts with. Both 
samples are fitted equally well by our code. A Gaussian fit to 
the Az/(1 + z) distribution shows an offset of 0.010 and scatter 
of cr — 0.028 for z < 1 . We get about 2% catastrophic outliers 
with |Az/(l + z)\ > 3cr - 0.084. This accuracy is a lower limit 
at bright magnitudes, while larger deviations are expected for 
fainter objects. We use a combination of the internally expected 
error and a lower limit that increases with magnitude and scales 
with (1 + z) as a more conservative error estimate. A z — 24 mag 
galaxy at z = 1.4 thus has a median error of cr z ~ 0.30. 

A cross-check shows our photo-zs to be in good agree- 
ment with the spectroscopic redshifts and the photo-zs from the 
COSMOS catalog (Mo basher eFaT1[2006l l below z = 0.4. Our 
code seems to underestimate redshifts at 0.4 < z < 1 .0, see upper 
panel of Fig.[T] At z > 0.75, the COSMOS catalog redshifts tend 
to be a little overestimated compared to the spectroscopic sam- 
ple, leading to a systematic difference of z ~ 0.08 ... 0. 1 with 
respect to our photo-zs. This trend seems to continue to higher 
redshifts (see lower panel of Fig.|TJ. 

Rest frame magnitudes are computed by convolving a red- 
shifted filter curve with the best-fitting template spectral energy 
distribution (SED) normalizing with the nearest observed band. 
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Figure 1. Upper panel: Comparison of spectroscopic and pho- 
tometric redshifts for HIROCS (filled cirles) and the COSMOS 
catalog (open squares). Lower panel: Comparison of HIROCS 
and COSMOS photometric redshifts up to z — 2. 
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Figure 2. //-band coverage of the COSMOS field and overdense 
regions in the redshift range 1.2 < z < 1.6. In the dashed region 
there is no //-band data at all, in the light grey regions it falls 
short of our target depth of H = 21 .4 (3cr). The IDs refer to 
Tabled 



To find galaxy clusters we use a method based on detect- 
ing excess density in the 3D galaxy distribution (Roser et al., 
in preparation) by computing the local density for all objects. 
For each object, the fractions of the photometric redshift prob- 
ability distributions p(z) of all its neighbours within 300 kpc in 
projected physical separation and a Az motivated by the average 
accuracy of our photometric redshifts are added. The resulting 
densities turn out to be rather insensitive to the exact value of Az. 
Only objects with 3<x overdensity compared to the average field 
are used to select cluster candidates. All structures with more 
than 6 overdense members within an aperture of 2' are picked 
as cluster candidates from a two-dimensional plot of overdense 
galaxies in redshift slices to avoid projection effects. Figure|2] 
shows such a plot for the redshift range 1 .2 < z < 1 .6 with the 
cluster candidates highlighted. Monte-Carlo simulations with 
COMBO- 17 data have shown this 3cr overdensity limit to be 
rather conservative (Falter et al., in preparation). Comparison 
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Table 1. Correlation of six z > 0.9 X-ray selected clusters from 
Table 1 of |Finoguenov et aL"1 (2006 ). ID refers to this table, zxmm 
is the redshift given there, Zhtrocs is the photometric redshift 
found by HIROCS. Msoo.x is the mass estimate based on X-ray 
luminosity, My itot the total rest-frame V-band luminosity. Mm/l,v 
is the mass estimate based on Mf iW . Masses are given in units of 
1O 13 M . 



ID 


Zxmm 


ZHIROCS 


Mjoo,x 


Af v>tot 


Mmil.v 


32 


0.90 


0.79 


18.64 ±0.52 


-25.0 


12.0 


66 


0.95 




1.81 ±0.45 






72 


0.90 


0.84 


5.34 ± 0.47 


-23.6 


3.2 


102 


1.10 


0.93 


2.60 ± 0.38 


-24.5 


7.4 


126 


1.00 


0.86 


6.87 ± 0.69 


-24.3 


6.1 


133 


1.15 




5.68 ± 0.64 







with the density maps of |Scoville et al. ( 2006[ l at z < 1.1 shows 
our overdense objects to trace all density peaks. A correlation of 
our HIROCS results with X-ray selected clusters at z >0.9 from 
Finoguenov et al. ( 2006 1 in our selection area shows that four 
out of six clusters are in common (see Table[T]i. The difference 
in redshift results from the different photo-z codes, see Fig.[T] 
For X-ray cluster 66 we find no counterpart in our //-selected 
object list. While this is a low-mass structure, X-ray cluster 133 
is also not detected by our algorithm since it has no members 
with a local overdensity > 3<x. Only few galaxies at its position 
are visible even if we lower the overdensity cut to lcr, despite its 
mass being above the mass of cluster 102 which we detect. All 
these X-ray selected clusters have redshifts Zhirocs < 1-2 and 
are thus not part of our sample. 



4. Results and discussion 

The first HIROCS cluster candidate sample for the COSMOS 
field comprises 15 galaxy cluster candidates with redshifts 

I. 22 < z < 1.55. They are listed in Table[2]with their positions 
and properties. Their significance with respect to the mean num- 
ber density of galaxies in this redshift slice is between 2.4<x and 

II. 4cr. Candidates #1 and #11 are removed from the final list 
due to their low significance, candidate #7 because it does not 
show a clear redshift peak and colour-magnitude plots show it to 
be most likely a chance projection of two structures. This leaves 
12 candidates in the final list. The number of cluster members 
is limited by the depth of the //-band data. Figure[3] shows an 
//-band image of HIROCS 095954.2+022924 (#13) at z = 1.45 
as an example. Most member galaxies are clustered in a region 
of about 1 □' on the sky. 

Some representative cluster members' SEDs from the tem- 
plate fitting process are shown in Fig. [4] Note how the z-band 
flux drops when moving up in redshift from z - 1.2 to z = 1-4 
(center column), placing a secure lower limit on the photomet- 
ric redshift. While we observe the very red, passively evolving 
galaxies expected to be present in galaxy clusters, there is a re- 
markably high fraction of galaxies with blue SEDs. This is also 
evident in colour-magnitude diagrams of the cluster candidates. 
Figure|5]shows restframe M280 - My colour-magnitude diagrams 
of all galaxies with H < 21.4 in the redshift slices [0.6, 0.8] 
(left) and [1.2, 1.6] (right). Overplotted in red are the members 
of a rich large-scale structure at a redshift of z ~ 0.7 (Guzzo et al. 



2007) (left panel) and members of the cluster candidates given in 
Table|2] (right panel). The strength of the colour bimodality de- 
creases with redshift, as observed for example by Cirasuolo et al. 
(2007 ). The fraction of red galaxies is considerably lower in our 
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Figure 3. //-band image 

HIROCS 095954.2+022924 (#13) at z 
galaxies are marked with white circles. 
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Figure 5. The restframe M280 - My colour-magnitude diagram 
in the redshift slices [0.6, 0.8] (left) and [1.2, 1.6] (right). Field 
galaxies are indicated as contours. Overplotted in red are the 
members of a known large-scale structure at z ~ 0.7 (left) and 
the HIROCS cluster candidates from Table[2] (right). 



high-redshift candidates than in the z = 0.7 structure. This could 
be due to a Butcher-Oemler type evolution, a colour-density ef- 
fect or a combination of both. 

On the HST ACS images many members exhibit disturbed 
morphologies and signs of interactions. This is seen for a lot of 
the blue and even for some of the reddest galaxies. Note that 
the term "red" only refers to colour and includes both passively 
evolving old stellar populations as well as dust-reddened objects. 

The total rest-frame V-band luminosity My to t of the candi- 
dates is in the range -25.8 < M Vtot < -24.2 with the bright- 
est cluster galaxies having rest frame luminosities in the range 
-24.3 < Mvbcg ^ -22.6. The brightest of these are consistent 
with present-day BCGs (Bahcall 2000) when applying a passive 
evolution correction factor based on the evolution of M* B from 



Bell 



lM/L,V 



et al.| (2004) with Zform = 3. Our mass estimates 
based on the total rest-frame luminosity are in good agreement 
with the X-ray based M500 masses, see Table[T] We assume mass- 
to-light ratios of 140 (z = 0.9) and 75 (z = 1.4) base d on the 
present-day mass-to-light ratio of 300 (Bahcall 2000) and the 
correction for evolution described above. Comparison with the 
mass estimates of the X-ray selected catalog shows that our algo- 
rithm detects structures down to at least M 5 oo,x = 2.6 ■ 10 13 Mq, 
the mass of X-ray selected cluster 102, at z = 0.95. Our lowest- 
mass z > 1.2 candidate HIROCS 095906.0+022726 (#15) has a 
total mass of M M j L V « 1.7- 10 13 Mq. At present we do not have a 
completeness limit estimate for our algorithm. Theoretical mod- 
els (ACDM, cr 8 = 0.8) predict a bout 12 clusters per .66 
at z > 1.1 with M > 5 ■ 10 13 M© ( |Bartelmann & White||2002 



M. Bartelmann priv. comm.). If we take M200 = 1-7 • M500 



( Finoguenov et al.|2 006) and assume M500 = Mm/l,v 11 of our 
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Table 2. First sample of HIROCS clusters in COSMOS with z > 1.2. An ID marked with "*" means that this candidate was rejected 
from the final list, see text for details. The cluster redshift z is the mean redshift of its members, cr, the error of the mean. N>3o- 
contains the number of members above the 3<x overdensity cutoff without correction for field galaxies, whereas N tot gives the total 
number of galaxies at the cluster redshift and position regardless of their local density with the average field density subtracted. 
They are selected to be within a circular aperture of 2x the median distance of the overdense members to the cluster center and to 
be in the redshift slice z c iuster ± 0.15. The significance of N tot with respect to the mean field density is given in "significance". M^bcg 
gives the total restframe V-band luminosity of the brightest cluster galaxy, My m the field corrected total luminosity in the restframe 
V-band down to the completeness limit. No correction for completeness is applied to My jto t- A^m/l, vgives a mass estimate based on 
Mvm, see text for details. Our 50% completeness limits are My = -21.4 at z — 1.3 and My = -21.8 at z — 1.5. 



ID 


Name 


KA L J 


DHL. L J 


z 


(T- 




XT 

Ntot 


!\jf 

'"V.BCG 


M vm si; 


*nificance 


A/I 

M M/LV 






(12000) 


(12000) 














[o-] 


[10 13 M o ] 


1* 


HIROCS 100035.6+020338 


150.148 


2.061 


1.22 


0.03 


7 


6 


-23.2 


-23.9 


2.4 


2.9 


2 


HIROCS 100220.5+020534 


150.586 


2.093 


1.23 


0.07 


9 


8 


-22.6 


-24.2 


4.5 


4.0 


3 


HIROCS 100212.2+014858 


150.551 


1.816 


1.24 


0.06 


9 


11 


-23.8 


-25.0 


5.3 


7.9 


4 


HIROCS 100223.7+014527 


150.599 


1.758 


1.24 


0.06 


10 


11 


-23.7 


-24.5 


3.5 


5.0 


5 


HIROCS 100054.5+024703 


150.227 


2.784 


1.26 


0.06 


7 


12 


-24.0 


-24.8 


9.0 


6.5 


6 


HIROCS 100221.8+021045 


150.591 


2.179 


1.31 


0.05 


8 


10 


-23.3 


-24.8 


6.4 


6.1 


7 * 


HIROCS 100023.9+024158 


150.100 


2.699 


1.32 


0.09 


17 


16 


-23.1 


-25.2 


9.1 


8.8 


8 


HIROCS 095957.9+024125 


149.991 


2.690 


1.33 


0.06 


8 


10 


-23.4 


-24.9 


7.3 


6.2 


9 


HIROCS 100103.2+024546 


150.263 


2.763 


1.40 


0.12 


6 


8 


-23.5 


-24.7 


6.6 


4.7 


10 


HIROCS 095931.8+023630 


149.882 


2.608 


1.41 


0.04 


6 


8 


-23.1 


-24.7 


4.4 


4.7 


11* 


HIROCS 095922.3+020654 


149.843 


2.115 


1.43 


0.04 


10 


19 


-23.4 


-25.8 


2.5 


13.1 


12 


HIROCS 095950.3+022014 


149.960 


2.337 


1.44 


0.07 


15 


19 


-24.3 


-25.9 


8.4 


14.3 


13 


HIROCS 095954.2+022924 


149.976 


2.490 


1.45 


0.05 


8 


10 


-23.6 


-24.9 


8.0 


5.2 


14 


HIROCS 095938.4+021939 


149.910 


2.328 


1.52 


0.06 


14 


18 


-23.0 


-25.3 


7.7 


6.8 


15 


HIROCS 095906.0+022726 


149.775 


2.457 


1.55 


0.05 


8 


8 


-23.2 


-23.9 


11.4 


1.7 




# 374590 
z = 22.98™° 
z= 1.51 +/-0.05 
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Figure 4. Representative SEDs of cluster candidate members in two redshift slices, z ~ 1.25 in the upper and z ~ 1.5 in the lower 
row. The left column shows very red SEDs whereas the middle and right columns show blue SEDs, indicative of recent or ongoing 
star formation. The redshift errors quoted are the results of the template fitting algorithm. 



12 candidates fall into this mass range, matching the predic- 
tions. Note that no contamination or completeness corrections 
have been applied. The cluster candidates we find are compa- 
rable in redshift to the most distant spectroscopically confirmed 
(McCarth y et al.|2007] > and photometrically selected (Cast e llano] 
|et al.|2006| > clustersknown today. These results indicate that with 
the full HIROCS galaxy cluster sample implications on cosmo- 
logical parameters and galaxy evolution will become assessible. 
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